Transport phenomena in three-dimensional system close to the magnetic quantum 
critical point: The conserving approximation with the current vertex corrections 



o 
o 



in 

(N 



c3 



C 

O 

i o i 

> 

O 

(N 

in 
o 

03 



o 
o 



X 



Sciichiro Onari, Hiroshi Kontani 1 , and Yukio Tanaka 
Department of Applied Physics, Nagoya University, Chikusa, Nagoya 464-8603, Japan, and 
1 Department of Physics, Nagoya University, Chikusa, Nagoya 464-8603, Japan. 

(Dated: February 2, 2008) 

It is known that various transport coefficients strongly deviate from conventional Fermi-liquid 
behaviors in many electron systems which are close to antiferromagnetic (AF) quantum critical 
points (QCP). For example, Hall coefficients and Nernst coefficients in three-dimensional heavy 
fermion CeCoIns and CeCus-zAu^ increase strikingly at low temperatures, whose overall behaviors 
are similar to those in high-T c cuprates. These temperature dependences are too strong to explain 
in terms of the relaxation time approximation. To elucidate the origin of these anomalous transport 
phenomena in three-dimensional systems, we study the current vertex corrections (CVC) based on 
the fluctuation exchange (FLEX) approximation, and find out decisive role of the CVC. The main 
finding of the present paper is that the Hall coefficient and the Nernst coefficient strongly increase 
thanks to the CVC in the vicinity of the AF QCP, irrespective of dimensionality. We also study the 
relaxation time of quasi-particles, and find that "hot points" and "cold lines" are formed in general 
three-dimensional systems due to strong AF fluctuations. 

PACS numbers: 72.10.Bg,71.27.+a 



I. INTRODUCTION 

Strongly correlated electron systems close to the quan- 
tum critical point (QCP) have stimulated much interest. 
Especially, heavy fermion compound CeCoInsji which is 
a three-dimensional metal close to the antiferromagnetic 
(AF) QCP, attracts much attention because of the non- 
Fermi-liquid normal state and the d-wave superconduc- 
tivity (T c = 2.3K). Moreover, recent experimental ef- 
forts have revealed the existence of interesting anoma- 
lous transport phenomena characteristic of the AF QCP. 
In the normal state of CeCoIns, for example, it is ob- 
served that the resistivity^ p oc T, the Hall coefficient 
Rh oc 1/T, and the Nernst coefficient v oc 1/T be- 
low 20K till T c = 2.3K. These behaviors are quite differ- 
ent from the normal Fermi-liquid behaviors, p oc T 2 and 
Rh oc T°. 

In Cell5, the maximum values of Rh and v are quite 
huge compared to values for high temperatures. The 
maximum value of Rh becomes about 30 times larger 
than the values at high temperatures. The value of v 
reaches — 1.0/iV/KT at 5K, whose magnitude is about 
1000 times larger than the values in usual metals. The 
temperature dependence of Rh is similar to the two- 
dimensional high-T c cuprates above the pseudo-gap tem- 
peratures, and that of v is similar to electron doped high- 
Tc cuprates, irrespective of signs. Moreover, the magni- 
tude of Rh and v in CeCoIns is much larger than that 
in high-r c cuprates. Similar drastic increases of Rh are 
also observed in CeCug-zAu^ and YbRh^Si^, which 
are three-dimensional heavy fermion compounds close to 
the AF QCP. 

The relaxation time approximation (RTA)£*£ has been 
used frequently in the study of transport phenomena, 
although its reliability for strongly correlated systems is 
not assured. According to the spin fluctuation theory, the 



relaxation time becomes strongly anisotropic The 
spots on the Fermi surface where takes the maximum 
(minimum) value is denoted by the cold (hot) spots in 
literatures. The ratio of the relaxation time at cold spots 
and hot spots r = r co id/Th t and the weight of the cold 
spots play an important role in the transport phenomena. 
However, in terms of the RTA, an unrealistic huge r (say 
r = 100-1000) is required to reproduce the experimental 
enhancement of Rh in CeCoIn 5 i£ If we assume that Rh 
is enhanced by this mechanism, Rh should be suppressed 
quite sensitively by a very small amount of impurity. In 
addition, when r > 1, the magnetresistance should be 
too large to explain the modified Kohler's rule, Ap/p oc 
R H /p 2 , which is observed in high-T c cuprates&iS and in 
CeCoInt- 

These anomalous transport phenomena close to the 
AF QCP are well reproduced by taking the current 
vertex corrections (CVC) into account. Actually the 
CVC is necessary to satisfy conservation laws. In the 
Fermi liquid theory, the CVC corresponds to the back- 
flow, which naturally arises from electron-electron corre- 
lations. Then, the CVC is indispensable to calculate the 
transport coefficients in the strongly correlated electron 
systems, where electron-electron correlations are dom- 
inant. For example, the modified Kohler's rule is ex- 
plained due to the CVC caused by the AF fluctuations^ 
The negative Rh in electron doped high-T c cuprates, 
which cannot be explained by the RTA because the Fermi 
surface is hole-like everywhere, is explained if we take the 
CVC. Moreover, it is not easy to explain the enhancement 
of v by the RTA because the Sondheimer cancellation^ 
makes v small. The enhancement of v in electron-doped 
(hole-doped) high-T c cuprates are caused by the CVC 
due to the AF (AF and superconducting) fluctuations. 

Until now, various non-Fermi-liquid behaviors of high- 
Tc cuprates have been explained by the spin fluctu- 
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ation model, such as the self-consistent renormaliza- 
tion (SCR) theorjiiii^ii^ and the fluctuation exchange 
(FLEX) theorj«iS*iii£. For example, an appropriate be- 
havior of the spin susceptibility and the AF correlation 
length £ (£ 2 oc 1/T) are obtained. Moreover, spin- 
fluctuation theories&i^ derive the relation p oc T 2 £ 2 , 
which is consistent with the non-Fermi-liquid behaviors 
of high-T c cuprates close to the AF QCP. Based on the 
spin fluctuation theory, Kontani et idJ-LLLtiidL have de- 
veloped a theory of transport phenomena, by focusing 
on crucial role of the CVC. This framework naturally 
reproduces the temperature dependence of transport co- 
efficients for high-T c cuprates and other 2D systems^ 
close to the QCP. 

In high-Tc cuprates^, the CVC plays an important 
role on transport phenomena. However, it is highly 
non-trivial whether the CVC is significant in three- 
dimensional systems, since the CVC totally vanishes in 
the dynamical mean field theory (DMFT) where d = oo 
limit is taken, 2 ^ The reason is that the irreducible four- 
point vertex T 1 in the DMFT becomes a local function, 
which cannot contribute to the CVC. Although it is gen- 
erally believed that the DMFT works well in various 3D 
systems with strong correlation^ the momentum de- 
pendences of the self-energy and the vertex corrections 
are significant near the QCP. However, we must con- 
sider the CVC which cannot be taken into account within 
the DMFT. To elucidate this issue, we study the role of 
the CVC based on the AF fluctuation theory in three- 
dimensional systems. 

The purpose of this paper is to examine whether non- 
Fermi-liquid behaviors in 3D systems can be explained 
by taking account of the CVC in terms of the FLEX 
approximation. In 3D systems, massive calculation re- 
source and time arc needed to perform the calculation. 
We show that striking increase of Rh and v can be ob- 
tained at low temperatures even in 3D by virtue of the 
CVC, which is consistent with experiments in CeCoIns 
and CeCu6-xAu x . The present study is the first micro- 
scopic calculation for the Hall coefficient and the Nernst 
coefficient in 3D with the CVC. We also find that in 3D 
systems, the hot and cold spots form point-like ("hot 
points") and line-like shape ("cold lines"), respectively. 
The CVC on the cold lines plays a major role for increas- 
ing Rh and v. 



II. FORMULATION 



A. Model 

We first introduce the three-dimensional Hubbard 
model, 

ft = - X X li ^ia c i° + 2 X U n ™ n i-° C 1 ) 

on a stacked square lattice with the Coulomb repulsion 
U and the intralayer hopping t\, t2, £3 and the interlayer 



hopping t z depicted in Fig. ^ Hereafter, we take t% = 1 
as a unit of energy. 
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FIG. 1: The lattice structure with the nearest intralayer hop- 
ping ti, the second one £2, the third one tz and the interlayer 
hopping t z 



Then, the dispersion is derived as, 

el = — 2ti(cos(fc x ) + cos(ky)) — 4i2 cos(k x ) cos(k y ) 

-2< 3 (cos(2fc x ) + cos(2fc y )) - 2t z cos(fc z ). (2) 

We apply the FLEX approximatio n 1 ^i 1 Ti 2 ^i 2 ^ where the 
Green's function, the self-energy and the susceptibility 
are obtained self-consistently. The FLEX approximation 
belongs to "conserving approximations" formulated by 
Baym and KadanofSZ^. The spin (xq) an d charge (x q ) 
susceptibilities are 



and the irreducible susceptibility Xa is 



(3) 
(4) 



T 



X° q (ui) = -jj X Gq+fchv, + e n )G fe (e„), (5) 



k . n 



where Matsubara frequencies are denoted by e n = (2n - 
l)7rT and oji = 21itT, respectively. 
The self-energy is given by 



E fc (e„) = — X G fc-<?( e « - vi)V q (u>i), (6) 

q,l 



where the effective interaction V q is 

V q (wi) = u 2 M + \xq("i) - x°K)) + u. (7) 

We calculate Green's function self-consistently with 
Dyson's equation, 



G fe (e„) 1 = ie n + /i - e° k - S fe (e„). 



(8) 
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The irreducible particle-hole vertex T 1 ,, , which satisfies 
r fefc'( e "' e »') = s ^k(t n )/SG k '(e n ,) is given by 

r 1 ^/ (e„,e„0 = V k _ k >{e n - e n /) (9) 

where the Maki-Thompson term is taken into account 
and the Aslamazov-Larkin term is omitted because the 
latter is negligible for the CVCpi^ 

In this paper we take N = N x x N y x N z = 64 x 64 x 32 
fe-point meshes and the Matsubara frequencies e n takes 
the value from -(2N C - 1)ttT to (2N C - 1)ttT with N c = 
256. 



B. conductivity 

In order to derive the transport coefficient, we begin 
with Kubo formula, 



kk era' 



l ™ K k*.k'*> + iS) 



, (10) 



where K(oj + iS) is analytic continuation of below K(iu>) 



dre^ 



TV{cL(T)c feCT (r) C + feV ,c feV ,})(U) 



and 



o = £^fc 

dk' 



EliashbergSS derived the conductivity in this way. By 
generalizing Eliashberg's theorjiiSiS the conductivity is 
obtained as 



Re[G*( £ )> fc:c ( e ) 2 } 



« fc (e) = V[e£ + ReE fc (c)], 



(13) 



where v k is the quasi-particle velocity (without Z- 
factor) and the retarded Green's function G k (e) is de- 
rived by the analytic continuation. The total current J k 
is given by the Bethe-Salpeter equation 

/°° df' 
— 7? 2 (fee,feV)|G^(e')| 2 J fe '(e'), 

which is based on the Ward identity. The irreducible 
four-point vertex T^ n (ke,k'e') is defined in Ref. |2^ Ac- 
cording to eq. ©, it is given by 



7^ 2 (A;e, fcV) = ^cotanh 



e' - e , e' 

tanh — 

2T 2T 



x2HmV k >_ k (e' -e + i5). (16) 



C. Hall coefficient 

When a weak magnetic field B is induced along z axis, 
the Hall coefficient Rh is given by 



Rh 



(17) 



The Hall conductivity a xy is obtained byi2i211 



A k (e) 



2tt V de 

Ac " x 

|ImG*(e)||G*( e )| 2 A fc ( ei . 
Wfex(e) 
+(x <-> y), 



(18) 



J fc ,( e )%^-^( e )^ (e) 



(19) 



where /(e) = l/(exp(e/T) + 1). 



D. Nernst effect 

Nernst coefficient under a weak magnetic field B along 
z axis and gradient of temperature along x axis is defined 
as 



Bd x T 



(20) 



(12) According to the linear response theoryai, the response 



function L 2 ^ is defined as 



T 



Jo 



dre^(T T jQ(r)j u (r = 0)). (21) 



The electron current operator j and the heat current op- 
erator j Q are given by 



k.a 



(22) 



(14) and 



j Q (lui) = lim V / dre^ T 

I (JL - A 

2\8t dr 
'Y ^ 



-J c kJ T ) c kA r ') 



(23) 



respectively. 

After the analytic continuation for eq. (|21[) . L xx and 



L X y are obtained as 32 



x {\Gl(e)\ 2 J kx (e) -Be[G%(e) 2 ]v kx (e)}&5) 
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L 2 X y(+i8) 



7fe(e) 



Q fc (0 



x|ImG£( £ )||G£( e )|X(e) 
-ImE fc (e), 



7fe(e) |Qfc(<0 x [«fc( £ ) x V 
fe' J 

x7? 2 (fc e ,fc'e')|G£,(e')| 2 Q v (e') 



7fe(«) 



(26) 
(27) 

(28) 



(29) 



where q k is the quasi-particle heat velocity <7fc(e) 
evfc(e). 

Using above expressions, v can be rewritten as 



-21 



Sa 



(30) 



(31) 



BT 2 a xx Ba xx 

where the thermopower S is given by 

1 T 21 
T 2 a ' 

■*- u XX 

III. RESULT 



Here, we show numerical results obtained by the CVC- 
FLEX approximation. We use filling n = 0.9 (n = 1 
corresponds to half filling) and the intralayer hopping 
parameters t\ = 1, <2 = — 1/6 and £3 = f/5 which repro- 
duce the Fermi surface of 2D high-X^ cuprate YBCO, 
and introduce the interlayer hopping t z which makes 
the Fermi surface three-dimensional. The Stoner factor 
as — max { J/Xq (a- 1 = 0)} represents the "distance" from 
the AF order (as = 1 corresponds to the boundary of 
the AF or the spin density wave (SDW) order) since the 
denominator of the static spin susceptibility is 1 — Ux°- 
We calculate for each t z with keeping as, by tuning the 
value of U as shown by solid line for as = 0.995 and 
dotted line for a s = 0.998 at T = 0.02 in Fig. El The 
"distance" from the AF order is considered to be same 
along these lines. Note that as < 1 is always satisfied in 
2D (t z = 0) at finite temperatures reflecting the theory 
of Mermin- Wagner 

Fig. Elshows the Fermi surface for t z = (2D), t z = 0.4 
(quasi-3D; q3D) and t z — 0.8 (3D). We see that three 
dimensionality becomes stronger as the value of t z in- 
creases. The fe-dependence of the static spin suscepti- 
bility xL ' s shown in Fig. 0] where the peak position 
is commensurate; (k x ,k y ) = (tv,tt) for t z — (2D), 
(k x ,k y ,k z ) — (7r,7r,7r) for t z — 0.4 (q3D) and incom- 
mensurate around (ir, it, it) for t z — 0.8 (3D). From these 
peak structures, we ensure that the AF fluctuations are 
dominant in these systems. In this case, the peak values 
of x s increase with t z as seen in Fig. [SJ This means that 
the present system approaches to the AF instability as 
the dimensionality changes from 2D to 3D. 



u 



^ n.,=«.99S 



tz 

FIG. 2: U against i z for as = 0.995 (solid line) and as 
0.998 (dotted line) at T = 0.02. 



To understand 3D structure of the relaxation time r, 
we show the fe-dependence of 7fe(= -ImE^fe = 0)) oc 
1/Tfc along the Fermi surface in Figs. 16181 where hot spots 
are depicted by circles and cold spots are illustrated by 
dotted circle. The bottom panels represent the momen- 
tum dependence of 7& along the Fermi surface for each 
k z . We see the hot spots exist on the plains of k z = 
and k z = ir/4 for t z = 0.4 and t z = 0.8, respectively. 

The volume fraction of hot spots decreases as t z in- 
creases, and for three-dimensional case (t z = 0.8) hot 
spots have point-like structure ("hot points"), which is 
consistent with the experimental result obtained by de 
Haas- van Alphen measurements on Celn3i2i Generally 
in 3D systems, nesting exists in only small parts of the 
Fermi surface, and hot spots form there. In general, 
the "hot lines'— where hot spots form line-like structure 
would not be appropriate in 3D systems. 

On the other hand, jk increases more sharply along fen 
direction than k z direction around minimum point of 7^ 
for t z = 0.4 and 0.8 as depicted in Fig. □ and Fig. 
In this sense, cold spots stretch strikingly along k z di- 
rection. Then, cold spots form line-like structure ("cold 
lines"). They are aligned perpendicular to the plain with 
hot spots. The formation of hot points and cold lines 
would be generally expected in three-dimensional sys- 
tems. To confirm the generality, we must study much 
more systems with various types of 3D Fermi surfaces. 

The calculated temperature dependence of the resistiv- 
ity p = l/a xx for t z = (2D), t z = 0.4 (q3D) and t z = 0.8 
(3D) with U = 9.4, U = 6.2 and U = 5.2, respectively 
(see Fig. are shown in Fig. [51 where the unit of p is 
hp - 2.4 x 10~ 6 fhTL for a bilaycr YBCO (lattice constant 
along c axis : a c — 5.8 x 10 -10 m) and 3.1 x 10 _6 f2m for 
CeCoIns (a c — 7.6 x 10 _10 m), respectively. The value of 
U is chosen to satisfy as — 0.995 at T = 0.02 for each 
case. In this case, the "distance" from the AF QCP is the 
same among three t z parameters. We see that, indepen- 
dently of the dimensionality, the resistivities with and 
without the CVC are proportional to the temperature. 
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FIG. 4: The spin susceptibility for t z = (left panel), 0.4 
(middle panel) with k z = 7T, 0.8 (right panel) with fc z = 7r , 
for q s = 0.995, T = 0.02. 



Y s 




FIG. 5: The temperature dependence of the peak value of 
spin susceptibility for t z — 0, 0.4, 0.8, in the condition of 
a s = 0.995 at T = 0.02. 



Then, the value of p is slightly enhanced by the CVC. It 
decreases as t x increases since the corresponding value of 
U is reduced. This T- liner behavior of the resistivity is 
consistent with experiments for the systems close to the 
AF QCP, such as two-dimensional high-T c cuprate and 
three-dimensional CeCoIn^. In detail, the resistivities 
with the CVC show sub-liner temperature dependence 
in low temperatures, which are also observed in the ex- 
perimental results for heavy fermion CeRhlns (private 
discussion). According to the SCR theoryfi^S the re- 
sistivity behaves as p cx T in 2D, and p oc T 3 / 2 in 3D. 
However, the SCR theory also predicts that p oc T for a 
wide range of temperatures even in 3D systems when the 
system is close to the AF QCP, which is consistent with 
the present numerical calculation. 

Next, we show the temperature dependence of the 
Hall coefficient Rh in Fig. ^3 where the unit of Rh 
is - 5.2 x lCT 10 m 3 /C for a bilayer YBCO (lattice 

constant along a and b axis : a a = = 3.8 x 10~ 10 m 
and that along c axis : a c = 5.8 x 10~ 10 m) and 
1.0 x lQ- 9 m 3 /C for CeCoIn 5 (a a = a b = 4.6 x 10- 10 m 
and a c = 7.6 x 10~ 10 m), respectively. In regard to 
the horizontal axis, unit of the temperature t\ is 4000K 
for YBCO. On the other hand, we estimate the nearest 
neighbor hopping t\ = 400K for CeCoIn 5 , because the 
experimetal data 2 shows that magnitude of Rjj begins 
to increase below 40K, which corresponds to O.lii in Fig. 
ITUl The Rh without the CVC, which corresponds to the 
RTA, is almost constant. However, independently of the 
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FIG. 7: 7 along the Fermi surface for t z = 0.4, as ~ 0.995 
at T = 0.02, where the hot spots and the "cold lines" are 
depicted by circles and dotted thick line, respectively. The 
trajectories of k\\ are depicted on the Fermi surface in k z — 
0, tv as examples. 



dimensionality, Rh with the CVC increases as tempera- 
ture decreases, which is consistent with experimental re- 
sults for high-T c cuprate^i and heavy fermion compounds 
(CeCoIn^, CeCug^Aa, 4 and YbRh 2 Si 2 5 ). Namely, the 
RTA cannot explain the strong temperature dependence 
of the Hall coefficient close to the AF QCP. Moreover, 
in 3D case hot spots take point-like shape, which means 
that the effective electron density for transport phenom- 
ena (n e ff) is large compared with two-dimensional case. 
Since RW A k is satisfied in the RTA, R c $ cannot 
become large in 3D systems. 

As a result, the CVC is indispensable to explain the 
behavior of Rh in 3D close to the AF QCP. In the 
present results, the maximum enhancement of Rh is 
given by R H (T = 0.02) /R H (t = 0.2) - 5 for t z = 0.8 
and as — 0.998. Rh should increase further if we cal- 
culate at lower temperatures. The shape of the Fermi 
surface in CeCoIns resembles to that of our model for 
t z = 0.8. However, to reproduce the experimental results 
in CeCoIns quantitatively, we have to study it based on 
the realistic band structures of CeCoIns Our calcu- 
lation shows that Rh is strongly enhanced by the CVC 
even in 3D systems, and its maximum value becomes as 
large as that in 2D. 

Here, we discuss the reason why Rh is strongly en- 



FIG. 8: 7 along the Fermi surface for t z = 0.8, as = 0.995, 
at T = 0.02 where the hot spots and the "cold lines" are 
depicted by circles and dotted thick line, respectively. 




°0 0.04 0,08 0.12 ' 0,16 0.2 

T 

FIG. 9: The resistivity p with the CVC (solid line) and with- 
out the CVC (dotted line) for t z = 0, 0.4, 0.8 as a function of 
T with a s = 0.995 at T = 0.02. 



hanced by the CVC in 3D systems based on the numeri- 
cal study. The general expression for a xy is given by^S^ 




where fcii is the component of k along the unit vector 
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ey(fe) = [e z x v k )/\v k \ which is in the fc^fcy-plane and 
parallel to the Fermi surface. \Jk\± = 



Jkx 



Jl y , and 



9 j is the angle between the total current J and the x 
axis. In this line integration fc-point moves anticlockwise 
along the Fermi surface around the k z axis. 

We see that for t z — 0.8, cold spots (where ~f k is small) 
form lines (cold lines) at the center of the Fermi surface 
along the z axis in Fig. |H| As shown in the last term 
of eq. i|32|) . main contribution for a xy is expected to 
come from the cold lines. We see that in Fig. El the 
momentum dependence of the absolute value of the total 
current \J k \ is quite similar to that of jk- 

In Fig. 12, we plot 9j along the trajectories of fc||. 
As references, we also plot the corresponding quantity 
without the CVC, 9 V for k z — 0-ir/2 and k z — Sti/A-tt as 
a thin line and a dotted thin line, respectively. As show in 
Fig. El Qv for k z = 0-n/2 and k z = Sir/A-ir decrease and 
increase, respectively, along the trajectories of ku . On the 
other hand, 9j has a non-monotonic change along these 
trajectories. Especially, for k z — 37r/4-7r 9j decreases 
contrary to the case of 9 V . We stress that the magnitude 
of d9j/dk\\ becomes larger than that without the CVC 
(<9#i;/<9fc||)around the cold lines. 

In Fig. 1181 we plot the momentum re- 



solved Hall conductivity a xy (k) 
= Z)fc cr xy(fc), where cr xy (k) 



defined by 
is given by 

^ y (fe) = -JSr ("§9 \lmGl(e)\\G^A k (e). 
The magnitude of a xy (k) takes large values around the 
cold lines , especially for k z = — 7r/2-7r/2. We should 
comment that for k z = 37r/4, ir the difference of the value 
of 9 V and 9j becomes n at the edge of the trajectory 
as shown in Fig. 1121 due to the CVC. At that time, the 
direction of J k is opposite to that of v k . In this case, 
strong AF spin fluctuation enhances the magnitude of 
72 2 (fee,feV) in the eq. (JTSJl for k - k' = ±(tt,tt,tt), 
±(— 7r, 7r, 7r), ±(7r, — 7r,7r) and ±(7r, 7r, — 71"). In this 
case, from eq. (|15|l . we can obtain J k ss v k + a k J k > 
(0 < a < l) 19 . This equation is easily solved as 



J k 



Vk + akVk' 
1 - a k a k > 



(33) 



In Fig. El we illustrate schematic behaviors of the 
quasi-particle velocity v k and the total current J k on 
the Fermi surfaces sliced at k z = tt (solid circle) and 
at k z = (dotted circle). We focus on the position 
of points A and B, which are connected by the nest- 
ing vector k — k = (— 7r, — 7r, tt). Here, we write the 
quasi-particle velocities at A and B as v k and v k >, re- 
spectively. We see that v k and v k > are antiparallel and 
\v k \ < \v k >\. Considering that a < 1, J k given in eq. 
(|33fl takes an opposite direction of v k . In the same way, 
the total current and the quasi-particle velocity at C are 
also antiparallel. This nontrivial behavior of J k has not 
been pointed out in previous studies for two-dimensional 
systems. This feature might induce an anomalous trans- 
port phenomenon. 



In Fig. 1151 we show the temperature dependence of 
R H for a s = 0.998 (solid line) and a s = 0.995 (dotted 
line) for t z = 0.4 and 0.8. We see that Rh with the 
CVC increases as as approaches unity. It seems that 
Rh tends to diverge as the system approaches the AF 
QCP. The reason can be understood by seeing Fig. 1161 
where d9j/dk\\ takes a large value at the cold spot which 
corresponds to the center of the fen axis. According to eq. 
(|32(l . this fact leads to the strong enhancement of a xy . 




"CM 



T 



FIG. 10: The Hall coefficient R H with the CVC (solid line) 
and without the CVC (dotted line) for t g = 0, 0.4, 0.8 as a 
function of T, with a s = 0.995 at T = 0.02. 



k?.=7C/4 




FIG. 11: Absolute value of the total current J with the CVC 
along the Fermi surface for t z — 0.8 and t z = (thin line) 
with as = 0.995 at T = 0.02. 



Finally, we discuss the Nernst coefficient v. It is known 
that v vanishes in a complete spherical system, which is 
called Sondheimer cancellation. 12,41 Although this can- 
cellation is not perfect in real anisotropic systems, the 
magnitude of v becomes small (~ InV/KT) in conven- 
tional metals. 

However, v is enhanced below T* (in the pseudo-gap 
region) for high-T e cuprates. The authors of Refs. |4l], 
l42l suggest a possibility that the vortex-like excitation 
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FIG. 12: the angle 0j of the total current J with the CVC 
along the Fermi surface for t z = 0.8, a s = 0.995 at T = 0.02. 
Thin line and dotted thin line approximately correspond to 
the angle 6 V of the current without the CVC for k z — 0-n/2 
and k z = 3ir/4-ir, respectively 




FIG. 13: u xy (k) along the Fermi surface for t z = 0.8, as 
0.995 at T = 0.02. 



emerges in under-doped high-T c cuprates to explain the 
enhancement of v in the pseudo-gap region. 

On the other hand, one of the present authosSi^ has 
shown that strong enhancement of v for high-T c cuprates 
is naturally derived based on the FLEX+T-matrix ap- 
proximation with the CVC. Furthermore, CeCoIn 5 also 
shows huge negative v below 20K, 3 which cannot be as- 
cribed to the vortex mechanism. Here, we aim to reveal 
the mechanism of the unconventional enhancement for v 
close to the QCP irrespective of the dimensionality. 

We estimate the renormalization factor (z) dependence 
of the transport coefficients, before showing the result of 
v. In the following, we will show that p and Rh are 
independent of z, and v oc z^ 1 in Hubbard model. Using 
the relation 3 ^ 

!/',„„ 1 f dS k de° k 1 



E = 

k 



(27T) 3 



dS k dk± — 



(2k) 



(27T) 



dS k de* k 
z\v k \ 
(34) 




FIG. 14: Schematic illustration to show that at A and B 
J k can direct opposite to v k , where the Fermi surface for 
k z — 0, 7r are depicted by dotted and solid circle, respectively, 
and J k and v k are by solid and dotted arrows, and the thick 
left-right arrows represent the coupling points. 
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FIG. 15: The Hall coefficient Rh with the CVC (upside) and 
without the CVC (downside) for t z — 0.4, 0.8 as a function of 
T with a s = 0.995 (solid line) ,0.998 (dotted line) at T = 0.02. 



where S k is the Fermi surface and k± represents the 
momentum perpendicular to the Fermi surface and also 
using the relation \G^(e)\ 2 ~ 7rz5(e-e* k )/j k (e* k = z(e° k + 
ReSfc(O) - n)) for 7 « T, we obtain 



dS k de* k ( df 



z\v k \ 



(2tt) 3 
cx (z)°. 

In the same way, a xy is given by, 

e 3 B f dS k del ( df 



--5- v kx J kx — (35) 

oe j c _ e , 7fe 



(36) 



4(2tt) 3 
oc (zf. 



z\v k \ 



dc 



z 

7 k 



A k ^ (37) 



(38) 
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FIG. 16: the angle 9j of the total current J with the CVC 
along the Fermi surface for t z = 0.8, k z = 7r/4, as = 0.995 
(solid line) and a s = 0.998 (dotted line) at T = 0.02. Thin 
line approximately corresponds to one without the CVC. 



We see that a xx and a xy are independent of z. Thus, we 
confirm that p = l/a xx oc (z) and Rh = H "I oc (z)° 

are independent of z. On the other hand, thermopower 
5 is given by 



S 



e 

{2ir) 3 Ta x 

en 2 k%T 
3(2n) 3 a x , 



f dS k de% ( df\ . 2 



2i 

de 



(4) : 



9 / "faJfa 



(39) 



(40) 

6 = 



z\v k \ dk± V I«fel7fe 



and using Q fc ~ q k 
Lll(+i5)/T 2 is given by 



ev k for T 



•"fcl7fc 

(41) 
(42) 
0, a xy = 



Be 2 



dSkde% ( df 



2(2tt)3Tj z\v k \ 
Be 2 f dS k de* k 



2(2tt) 3 T J z\v k 
Be 2 TT 2 k 2 B T f dS k d 



de 

df 
de 



4- ) eSUft(4)4 (43) 

<k 



6(2tt) 3 
-l 



- c k 

(4) 2 



4" 



z|« fc | <9fc x Vl«fcl7fe7 E=0 



where AJ£(e) is defined by 
^fc(e) = 7fc(e) j«fc(e) * [« fc (e) x V] a 



Jfc(e) 



7fe(e) 



(45) 
(46) 



(47) 



Thus, we obtain 



-,- ^ z . We must 

Baxx Ba xx 

consider z in detail to calculate v, because z is much 
smaller in heavy fermion systems. From the experiment 
of the de Haas-van Alphenii, we can estimate that the 
effective mass m* ~ lOOmo (mo is the bare electron mass) 



and the mass obtained by the band calculation mj, ~ 2irio 
in branches f3. Then, the mass enhancement factor z~ 1 
is given by z _1 = m* /nib ~ 50. 

Although the FLEX can describe various critical phe- 
nomena near the AF QCP, the mass enhancement is not 
completely explained with the FLEX in heavy fermion 
systems. The reason is that local correlations are not 
fully taken into account in the framework of the FLEX, 
because the vertex corrections in the self-energy are not 
included. According to Ref. 0, we separate the self- 
energy into the "local part" and the "non-local part" . 
In this case, total renormalization factor z is obtained 
as z = zqz*, where zq is the local renormalization factor 
which cannot be included in the FLEX and renormal- 
ization factor z* = (1 — (dT./dui))^ 1 is obtained by the 
FLEX. To fit the total renormalization factor to the ex- 
perimental results (z~ x = 50), we use zq = 3/50, because 
z* ~ 1/3 in our calculation. 

We show the obtained temperature dependence of the 
Nernst coefficient v in Fig. El where solid and dotted 
lines correspond to v with and without the CVC, respec- 
tively, for t z = 0, t z = 0.4 and t z = 0.8. In this figure, 
we chose the parameters for heavy fermion CeCoIns, i.e., 
a a = a b = 4.6 x 10~ 10 m, and then, fcB °° a " - 28nV/KT 
has been multiplied as a unit of calculated value, and the 
local mass enhancement factor z^ 1 has been also multi- 
plied. 

We see that v without the CVC is almost constant, 
and v with the CVC shows enormous increase at low 
temperatures, especially in strong three-dimensional case 
(t z — 0.8) where the Fermi surface is similar to that of 
CeCoIns. Then, the temperature dependence of v re- 
sembles to that of Rh . This temperature dependence of 
\v\ is consistent with the giant Nernst effect in CeCoIns 
(H - 1/iV/KT for T = 5K)£ Here, we discuss the rea- 
son why the magnitude of v becomes large. A k and A' k 
given in eqs. Ijl9(l and (|28J) are rewritten as 



A fe (e) 



d 



bfe(e)|± ( Jk{e) x gTJ- J fe( e ' 
\v k {e)\ ± \J k (e)\ 2 



d9j(k) 



0k 



(48) 
(49) 



A'k(e) = lk(e)\v k (e)U [ QM x I (50! 



\v k (e)l 



{Qkx^kx ~\~ QkyJky^) 



dk\\ 7fe(e) 
dOj{k) 



dk n 



+(Q fc x Jk)z JLi J\d]± 
oh V 7fc 



(51) 



Here, Q k is the total heat current with the CVC. We 
stress that Q k ~ ev k and Q k is not parallel to Jk 
when the AF fluctuations are strongiSl In this case, v 
is strongly enhanced due to the second term of A' k in 
high-T c cuprates4£ We expect that the same mechanism 
will give the enhancement of v in the present study for 
three-dimensional case. In NCCO, v is enhanced by the 
CVC due to the AF fluctuation below 300K, whereas 
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the increment of v for LSCO is brought by the CVC 
due to both AF and superconducting fluctuations below 
T* ~ 150Ki2i Because of the relation v oc 7~ 1 , v is 
proportional to p^ 1 for a fixed as- This fact would con- 
tribute to the enhancement of v in 3D shown in 
Figs. El El We show the temperature dependence of 
v for both a s = 0.995 and 0.998 in Fig. CHI v with the 
CVC increases as as approaches unity, which resembles 
to the behavior of Rh • Consequently, the magnitude of 
v increases almost divergently in the vicinity of the AF 
QCP. 



QiV/KT] 




FIG. 17: The Nernst coefficient v with the CVC (solid line) 
and without the CVC (dotted line) for t z = 0, 0.4, 0.8 as a 
function of T with a s = 0.995 at T = 5K 
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FIG. 18: The Nernst coefficient v with the CVC (upside) and 
without the CVC (downside) for t z = 0.4, 0.8 as a function 
of T with as = 0.995 (solid line) and 0.998 (dotted line) at 
T = 5K. 



IV. CONCLUSION 

We have calculated microscopically the resistivity p, 
the Hall coefficient Rh and Nernst coefficient v for three- 
dimensional Hubbard model close to the AF QCP based 
on the Fermi-liquid theory. This is a first microscopic 
calculation for the Hall coefficient and the Nernst co- 
efficient with the current vertex corrections (CVC) in 
the three-dimensional system. In two-dimensional sys- 
tems, it is established that the CVC plays an important 
role when the AF fluctuations are strong. On the other 
hand, the CVC vanishes completely in infinite dimension 
d = co. Thus, it is a very important theoretical issue 
to clarify whether the CVC is significant or not in three- 
dimensional systems. We find that the CVC influences 
crucially on various transport phenomena in both two 
and three-dimensional systems close to the AF QCP. 

We have shown that the magnitude of Rh and v is 
strongly enhanced with the decrease of temperature due 
to the CVC. These strong temperature dependences in 
the Hall coefficient and the Nernst coefficient come from 
the difference between the direction of the total current 
( Jfc) and that of around the cold spots. The difference 
of directions increases as the temperature decreases near 
the QCP, which can be expressed in terms of the effec- 
tive curvature dQj/dk\\ of Fermi surface obtained by the 
direction of Jk ■ The obtained values of Rh at the lowest 
temperature (T = 0.02) is more than 5 times larger than 
those at high temperatures for three-dimensional system 
(t z = 0.8). This result is qualitatively consistent with 
experimental results in various three-dimensional heavy 
fermion systems close to the AF QCP, such as CeCoIn5 
and CeCug_ a; Au 2 ;. This strong enhancement cannot be 
explained by the relaxation time approximation (RTA). 

In the present paper, we also studied the momen- 
tum dependence of relaxation time in three-dimensional 
systems due to strong AF spin fluctuations. In two- 
dimensional systems, it is known that hot spots and cold 
spots take line structures along z axis, as the systems 
approach to the AF QCP. In three-dimensional systems, 
we find that hot spots become point-like ("hot points") 
while the cold spots remain to take line structures ( "cold 
lines" ) . The emergence of hot points and cold lines is ex- 
pected to be general in three-dimensional systems close 
to the AF QCP. Transport phenomena are mainly deter- 
mined by the cold spots. The area of cold spots in the 
phase space plays an important role. We find that the 
CVC around cold spots produces strong enhancement of 
Rh and v. We emphasize that the strong enhancement 
of Rh and v comes from the effective curvature of the 
Fermi surface, dOj/dkn, enhanced by the CVC on cold 
lines, as shown in Fig. [S| Note that obtained results for 
Rh and v without the CVC, which corresponds to the 
RTA, are almost temperature independent. 

In future, we will perform a quantitative study for 
the transport phenomena in CeCoIns and CeRhlns, 
using a realistic band structure predicted by band 
calculations £IL2iL In the present paper, signs of Rh and 
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v are opposite to actual experimental results. We ex- 
pect that this discrepancy can be resolved by taking into 
account a proper band structure. 
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